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The reheating process in inflationary universe models is considered as an out-of-equilibrium mix- 
ture of two interacting and reacting fluids, and studied within the framework of causal, irreversible 
thermodynamics. The evolution of the temperature and the decay rate as determined by causal 
, thermodynamics are estimated at different stages of the process. A simple model is also used to find 

■ the perturbations of the expansion rate, including the possibility of damped oscillations. 



I. INTRODUCTION 

> 

The reheating period is a key ingredient in many inflationary models of the very early Universe. During reheating, 
most of the matter and radiation of the Universe is created via the decay of the inflaton field, while the temperature 
grows by many orders of magnitude. Traditionally, the idealized picture of a very quick decay in which the products 
reach equilibrium immediately is assumed JjJ. However, recent quantum field theoretical investigations g| - |Tc| l 
indicate that the reheating period is characterized by complicated nonequilibrium processes, the main characteristics 
of which are initial, violent particle production via parametric resonance ('preheating'), with a highly nonequilibrium 
' distribution of the produced particles, subsequently relaxing to an equilibrium state. 
Oh' In this paper we aim at obtaining a phenomenological understanding of the reheating process within a model of 
two interacting and reacting fluids. We regard this approach as complementary to the above-mentioned quantum 
field theoretical studies. The model involves a fluid with the equations of state for matter (modelling scalar field 
oscillations about the true ground state in the reheating) that decays into a relativistic fluid. The implications of an 
intermediate decay into other massive (bosonic) particles that does not explicitly occur in this model, are assumed 
to be describable with the help of a large, effective decay rate of the initial into the final component of the entire 
process, leaving the detailed microphysical study of this epoch to quantum field theoretical investigations which are 
beyond the scope of the present paper. 



While these simplifications may appear drastic, they open the possibility of studying the backreaction of the decay 
process as a whole on the entire cosmological dynamics, including the behaviour of the scale factor. We believe this 
to be a main advantage of our approach since most of the quantum field theoretical calculations do not even take into 
account the expansion of the universe. 

A two-fluid description of the expanding Universe is necessarily dissipative, even if the components are assumed to 
be intrinsically perfect fluids. This is also true for nonreacting fluids (conserved particle numbers), where the different 
cooling rates of the subsystems produce an effective, entropy generating bulk viscous pressure of the system as a whole 
{□J , Jl2[ . In fact conservation of the particle number is only a very special case, particularly at high energies. The 
processes we are interested in are characterized by particle decay and production. 
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As has been shown recently, any deviation from detailed balance in the decay and inverse decay reactions in the 
expanding Universe gives rise to additional bulk pressures Jl3| . In the present case of strong perturbations of the 
detailed balance (the reactions are predominantly one-directional with the inverse processes largely suppressed), one 
expects considerable bulk viscous pressures that characterize the deviations from equilibrium of the cosmological 
fluid as a whole. To describe this noncquilibrium process, we shall resort to the well-known Israel-Stewart theory of 
transport processes H] which, because of its causality and stability properties p5|, has been repeatedly applied in 



the cosmological context - see |16j], |I7| , |18 and references therein. 

In this theory, generalized 'fluxes', like a bulk pressure, become dynamical degrees of freedom on their own and 
have associated relaxation times. Our point of view here is to regard the cosmic substratum during the reheating 
period as a causal, dissipative fluid, relaxing to equilibrium. (By "causal" we mean that dissipative signals propagate 
only at subluminal speeds). 

Bearing in mind that the Israel-Stewart theory was derived for small deviations from equilibrium (for attempts 
to apply it to far-from-equilibrium situations see JL7| , Jl8[|), we point out that in the present context it is the high 
'preheating' particle production rate that provides a 'creation' contribution to the 'effective' bulk pressure which may 
be much larger than the conventional, thermodynamic, dissipative bulk pressure. While the latter is considered to 
be small in the applications of this paper, the 'effective' bulk pressure, that determines the reheating temperature, is 
not. 

Section II establishes the basic relationships concerning both particle and energy-momentum non-conservation in 
our two-fluid system, and presents an expression for the evolution of the temperature of the overall fluid. Section III 
introduces the causal transport relation for the thermodynamic, dissipative bulk stress, and deduces the corresponding 
equation for the evolution of the Hubble parameter. Section IVA solves the latter equation at three different stages 
of the reheating by assuming a very short relaxation time, and determines the temperature evolution in each of 
them. The qualitative description of the dynamics is determined essentially by a single thermodynamical quantity, 
the dissipative contribution to the speed of sound. If this is larger than the adiabatic contribution to the sound speed 
at the beginning of reheating, we show that the temperature rises rapidly to a maximum (reheating) temperature, 
which we estimate. This condition for rising temperature is equivalent to a growth in the total particle number 
density, which is reasonable in the initial stage of reheating, when coherent oscillations of the scalar field lead to 
huge production of particles. Thereafter, the temperature falls, but less rapidly than in the non-dissipative case. In 
Section IVB, we construct a simple model to calculate the perturbations of the expansion rate due to causal viscous 
and reactive effects. The model includes the possibility of damped oscillations in the beginning of reheating. Finally, 
section V summarizes our conclusions. 

Units have been chosen so that c = &b = 7i = 1. 



II. THE TWO-FLUID MODEL 



Let us assume the energy-momentum tensor T lk of the cosmic medium splits into two perfect fluid parts: 

rjiih rj-ilk rj-iik (1) 

with (A = 1, 2) 

T\ k =p A u i u k +p A h ik . (2) 

Pa is the equilibrium pressure of species A. For simplicity we assume that both components share the same 4- velocity 
it 1 , with projection tensor h lk = g %k + u % u k . The particle flow vector N A of species A is 

N\ = n A u l , (3) 

where n A is the particle number density. We are interested in situations where neither the particle numbers nor the 
energy-momentum of the components are separately conserved, i.e. particle interconversion and exchange of energy 
and momentum between the components are admitted. 
The balance laws for the particle numbers are 

N A;i = n A + ®n A = n A T A , (4) 

where = uf. is the fluid expansion and T A is the rate of change of the number of particles of species A. There 
is particle production for T A > and particle decay for T A < 0. For T A — we have separate particle number 
conservation. 
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Interactions between the fluid components amount to the mutual exchange of energy and momentum. Consequently, 
there will be no local energy-momentum conservation for the subsystems separately. Only the energy-momentum 
tensor of the system as a whole is conserved. Denoting the loss- and source-terms in the separate balances by t A , we 
may write 

n% = -A , (5) 

implying 

p A + 6 (p A + Pa) = Uif A , (6) 

and 

{pA+PA)u i +PA,kh ik = -h{t%. (7) 

All the considerations to follow will be independent of the specific structure of the t\. In general, there are no 
limitations on the strength or the structure of the interaction. 
Each component is governed by a separate Gibbs equation: 

T AdsA = d (£±) +PAd (±) . ( 8 ) 



n A ) \n A/ 

Using (|j) and @) one finds for the time evolution of the entropy per particle 

n A T A SA = Ui? A - (p A + pa) ^a • (9) 

With nonvanishing source terms in the balances for ua and pa, the rate of change of entropy per particle is nonzero 
in general. Below we shall deal with the special case that the terms on the right hand side of (Q) just cancel. 
According to (|^) the condition of energy-momentum conservation for the system as a whole, 

(Tf +T 2 k ), k =0 , (10) 

implies 

A = -A • (ii) 

There is no corresponding condition, however, for the particle number balance as a whole. Defining the integral 
particle number density n as 

n = ni + n 2 , (12) 

we have 

n + On = nT , (13) 

with 

nT = njTi + n 2 T 2 . (14) 

F is the rate by which the total particle number n changes. We do not require T to be zero, since total particle number 
conservation is only a very special case, especially at high energies. 

From now on we assume that the source terms on the right hand side of (^) cancel among themselves, i.e., that 
the entropy per particle of each of the components is preserved. The particles decay or come into being with a fixed 
entropy sa- This adiabaticity condition amounts to the assumption that the particles at any stage are amenable to a 
perfect fluid description. With sa = in @ one has 

Uit A = (p A + pa) ■ (15) 

This relationship establishes a link between the source terms in (Q) and (^) which originally are independent quantities. 
The simplifying assumption sa = takes into account the circumstance that the production process itself is the main 
source of entropy production, while dissipative processes within each of the separate components are less important. 
The cosmic fluid as a whole will be considered, however, as dissipative (see below). 
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Combining (11) and ( |l5| ) one has 

uit\ = (pi +px) Ti = -Uit\ = 
which provides us with a relation between the rates Ti and T 2 

Pi + Pi 



(pi +p 2 )r 2 , 



Pi + Pi 



Use of this relation in ( |l4| ) yields 



nT = niTihi 



1 

hi 



1 

~h 2 



(16) 



(17) 



(18) 



where Ha = (pa + Pa) In a arc the enthalpies per particle. Total particle number conservation, i.e. r = 0, is only 
possible if h\ — h 2 . 

As is shown elsewhere jl2| , jl3| , a system of two fluids, each of them perfect on its own, is dissipative in general 
and may be characterized by an energy-momentum tensor 

T lk = pu l u k + {p + n) h tk . (19) 

The equilibrium pressure p of the total system and the energy density p are assumed to obey equations of state 

p = p(n,m,T) (20) 

and 

p = p(n,m,T) , (21) 

where T is the equilibrium temperature of the system as a whole, defined by (cf. @, @) 

pi (rai,Ti) + p 2 (n 2 ,T 2 ) = p(n,ni,T) . (22) 

It is worth mentioning that there does not exist a corresponding relation for the pressures. The partial pressures of 
the components do not add up to the equilibrium pressure in general. The difference between the sum pi (m,Ti) + 
Pi {tii, T 2 ) and the equilibrium pressure p (n, T) contributes to the thermodynamical, dissipative bulk pressure n. A 
further source of it is deviations from detailed balance, i.e. contributions due to Fa 7^ p3|. 
The behavior of the equilibrium temperature T of the system as a whole is governed by]13| 



I = - ( e-r)**— *L 

T drp Tdrp 



where the abbreviations &rf — df/dT and 

pT-piTi -p 2 T 2 
Td T p 

Pi +Pi (n 2 n\ 
( — Pl 

V n n 




(23) 



Td T p 



Pi 



1 



1 



nxhx n 2 h 2 



(24) 



were used. The term (24) takes into account deviations from classical gas behavior. It vanishes for pa = tiaT, i.e. 
for a mixture of classical gases. For a mixture of a classical and a quantum gas, or for a mixture of fermions and 
bosons, it will be nonzero in general. We shall restrict ourselves to the case of two classical fluids, i.e. to the case 
(T/T)* = 0. This simplifying assumption is in line with the general restriction of our approach concerning the detailed 
microphysics, pointed out in the introduction. 

The temperature evolution equation ( |23| ) suggests that we define an 'effective' bulk pressure 



e 



p . 



(25) 



so that (E3J) in the classical case may be written as 
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: r s^i-^l. {26) 

1 o T p lo T p 

Our main concern in the following sections will be to study the influence of the different parts of the entropy producing 
'effective' bulk pressure n ef , on the cosmological evolution. 

The 7r c// -term in ( p6] ) describes the deviation from the adiabatic temperature behaviour. It is obvious that any 
TT eff < on the right-hand side of (|26| ) yields a positive, i.e., 'reheating' contribution to T/T, counteracting the first 
term on the right-hand side of ( p6| ) that simply describes the adiabatic coolin g du e to the expansion. 

As we shall discuss below the nonequilibrium term on the right-hand side of (|26|) may overcompensate the adiabatic 
term during the initial 'preheating' stage. 

The effective bulk pressure ([25]) consists of the conventional, thermodynamic, dissipative bulk pressure tt and a 
creation part — pT/Q. For T > G, a condition that is expected to be fulfilled during the initial stage of reheating, the 
creation part gives the dominant contribution to the effective bulk pressure. Even for a small or vanishing thermody- 
namic viscous pressure tt the entropy producing nonequilibrium parts on the right-hand side of the temperature law 
( pri| ) may overcompensate the adiabatic part, provided only that the production rate T is sufficiently high. In other 
words, even the production of particles with an equilibrium distribution, equivalent to the possibility of a perfect fluid 
description, gives rise to 'reheating' and entropy production. We shall consider this to be the dominant part of the 
entropy production during preheating. It is the particle production process itself which is connected with entropy 
production, simply through the enlargement of the phase space. There are additional entropy producing contributions 
due to the fact that in reality the particles will deviate from equilibrium. These contributions will be subject to a 
causal transport equation in the following section. 



III. THE CAUSAL EVOLUTION EQUATION 



In the Israel-Stewart second order theory of irreversible processes, the viscous pressure tt is a dynamical degree 
of freedom on its own 113]. Instead of applying the much more involved full Israel-Stewart theory we shall restrict 
ourselves to the so-called truncated version [O , [jl8| , [jl9| of this theory, since the latter already captures the essence 
of noninstantaneous relaxations. While the full and the truncated theories may disagree significantly if applied to 
far-from-equilibrium situations |l7j , |T§| ] , JlS| ] , they are expected to provide similar results near equilibrium p"8| ] . Now, 
reheating is anything but close to equilibrium. But as already mentioned, the main contribution to the entropy 
production during the decay process may be traced back to the second term in the expression ( |25| ) for the effective 
bulk pressure. The thermodynamic, dissipative bulk pressure tt may be regarded as a small perturbation under these 
circumstances and is supposed to obey the truncated causal transport equation 

TT + TTT = -(6 . (27) 

The quantity r denotes the relaxation time associated with the thermodynamic, disipative bulk pressure, i.e. the 
time the system would take to come to equilibrium (perfect fluid behaviour) if the generalized 'force' - in this case 
the expansion <d - were suddenly turned off. It may be related to the bulk viscosity £ by Jjl| (see also ITsf ) 

^=c 2 (p + p) 7 c 2 <l-c 2 s , (28) 

where c s is the adiabatic sound speed and cj, is the bulk viscous countribution to the dissipative speed of sound v, 
given by v 2 = c 2 + c 2 . Thus 

TT + TTT = —BclpjT . (29) 

Here and throughout we will use the abbreviation 

P 

where 7 is not assumed constant. 

The Einstein field equations for a spatially-flat Robertson- Walker spacetime are 

3H 2 = up , (30) 

where k is Einstein's gravitational constant, and 
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H 



±K(p + p + 7r) , 



(31) 



where H = ^8 = R/R and R is the cosmic scale factor. By eqs.([30]) and (|3l|), the thermodynamic bulk pressure may 
be written as 



ktt = -3jH 2 - 2H 



To calculate tt one first has to determine p via 



dp. . dp- 
P= d^ n+ df T 



Using the balance law ( |13| ) for n, and the evolution equation (|23|) for T, we find 

p = - (3if - T) (p + p) c 2 - 3Hir 
with the adiabatic sound velocity c s given by 



fdrp' 
\d T p) 



dp 

dp 



n dp T {pTp) 



isentropic 



(p + p) dn (p + p) d T p 



(32) 



(33) 



(34) 



(35) 



By differentiating (^Tj), we get 



2 H H ( d TP \ 

+ (3H-T) (p + P )c 2 s . 
Using ([32]) and (p^) in (p9|), the evolution equation for H becomes 



3H (p + p) 



drP 
8tp 



(36) 



& f^ dT P\^ 9 rr ( dT P 



H 3 

i^ + 2 7 = °- 



(37) 



Equation ( p7| ) is the causal dynamical equation in general form. We now specify the equations of state implicit in 
p and p. We assume that fluid 1 is described by the equations of state for nonrelativistic matter, i.e. 



while fluid 2 is relativistic: 



Consequently (B7I) reduces to 



pi = nim + \n\T , p\ — n{T , m 3> T , 



p 2 = 'in 2 T , p 2 = n 2 T . 



(38) 



(39) 



H 


( 5nx - 


f 8n 2 ^ 


f + 3 { 




V ni -{ 


- 2n 2 J 





{ni + n 2 ) ri\m 



(5rti + 8n 2 )nim 



(ni + 2n 2 ) (riim + 3n 2 T) 2 
H 3(n im + 4n 2 T) 



l4\ H 



8(ni + 2n 2 )(n\m + 3n 2 T) 



Q 



H 2 2(mm + in 2 T) 



(40) 



where Q = |Ti | > was assumed, i.e. the nonrelativistic component decays. This is the causal evolution equation to 
be solved. 
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IV. CAUSAL DYNAMICS OF REHEATING 



The equation ( [40| ) is a complicated nonlinear second order equation. To solve it, one could resort to numerical 
integration. However, this is strongly dependent on initial conditions and on the form of the decay rate Q, and it 
cannot readily give an idea of the overall dynamical features implied by the equation. Furthermore, part of the purpose 
of our thermodynamical approach is to avoid detailed complexities and to aim for an overall qualitative understanding 
arising from the constraint of causality. 

In Section IV. A, we use eq.([40|) to estimate the temperature and decay rate during reheating, neglecting the small 
perturbation of the expansion rate due to thermodynamic viscous effects. (Notice that the gravitational field equations 
'feel' only the thermodynamic viscous pressure.) It turns out that under this condition the reheating dynamics may 
be discussed in terms of the ratio Cb/c s of the dissipative to the adiabatic contributions to the speed of sound. Then 
in Section IV. B, we calculate the perturbation of the expansion rate in a simple case. 



A. Temperature and decay rate 



It is reasonable to assume that the expansion is approximately governed by non-viscous effects, and that the latter 
can be treated as a back-reaction. In this approximation, we neglect the terms multiplying r in (40), and arrive at 
the equation governing the expansion rate: 



H 3(nim + 4n 2 T) 
TP + 2(niTO + 3n 2 T) 







(41) 



We solve this for three different stages of the decay. Then we use ( ^0| ) to determine the corresponding decay rate Q. 
The evolution of the temperature is given by 



T 



T 6(ni+2n 2 )T 



-mr(pi 



P2 



nimQ] 



(42) 



as follows from (Eq) and from 



An T 



which is a consequence of ( |1S| ) and the equations of state fl38| ) and (39). In the perfect fluid limit n = Q = one 
recovers T/T = — 2H for n\ 3> U2, i.e. the nonrelativistic case, while for radiation [n\ <C n-i) the well-known 
behaviour T/T = —H is reproduced. Assuming that \ir\ <C p\ + P2, in agreement with the applicability conditions of 
the Isreal-Stewart theory, equation (E2h may be written as 



T 

- ps -2H 
T vni 



m + n 2 



2n 2 



1 



r 

3H 



(43) 



It follows that the temperature rises when T > 3H. By (|i3|), this is equivalent to h > 0, i.e. to a growth in the total 
particle number density (cf. j2(J). As discussed in Section I, such abundant net particle creation is expected to occur 
in the initial nonrelativistic stage of reheating. 



(i) Nonrelativistic regime 

In this regime the fluid is dominated by massive particles, so that m ^> 712, J ~ 1. Then ([y]) reduces to 



H+%H 2 



=> H re — 
3t 



(44) 



Inserting this solution back into (HQ), we find that the decay rate at the beginning of reheating is 



Q « f c 2 b H 



Cb 



H where c„ 



5 T 

3 TO 



(45) 
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It is a specific feature of the present first-order approximation that it fixes the rates Q and T which, in our general 
setting, are input parameters. Q and T are positive and proportional to H. Moreover, our approximation relates 
the bulk viscous contribution to the sound speed, to Q and T. The exact value of V is fixed by the ratio c^/cj.. Since, 



according to fl45|), Cj«l,a value « 1 is admitted by (|28|). This amounts to T ^> H, i.e., violent particle production 
as required for 'preheating'. A sufficiently high value of T also determines the entropy production. Neglecting second- 
order terms in 7r, the entropy flow vector is given by S a ~ nsu a , where s is the entropy per particle. For large T, 
neglecting the change in the entropy per particle, the entropy production density is approximately 



S"„ 



3(S 



sH . 



Equation ( f42| ) implies that at the beginning of the reheating, the rate of temperature change is given by 

2" 



T 



H . 



(46) 



(47) 



In the initial stage, a very large rate of creation of particles can lead to a growth in the net number density and thus 
in the temperature. By (47), this is again equivalent to the effective dissipative contribution q, to the sound speed 
exceeding the adiabatic contribution c s : 



n > T > <^ c b > c s 



(48) 



Equation ( [47| ) then implies that initially the temperature rises extremely rapidly (see also |20| ] ) . It reaches a maximum, 
which we could call the reheating temperature T re h, following standard terminology, and then decreases as Cb falls 
below c s . By (45), the reheating temperature is given by 



T, 



reh 



hi 171 



(49) 



It follows from ( fig ) that there is high entropy production in this regime. Alternatively, this regime may be char- 
acterized in terms of the effective bulk pressure ( p5[ 
to 



For T > SH the condition |7r| p± + p 2 — p leads generally 

r 



P 



Using the relations (^5|) one gets 



equivalent to 



p 



3H 



-—nmct 
5 



(50) 



(51) 



T, 



reh 



p 



For large particle production the condition |7T e ,, | 3> p is fulfilled. The reheating temperature 
as the temperature for which |7r e/ , | ~ p ~ p±. 



(52) 

may be understood 



(ii) Intermediate regime 



Here the energy densities of both components are comparable, i.e. n\m w 3n 2 T so that 7 




and consequently 



equation (|iT|) simplifies to 

H+lH^Q => H^^. (53) 

Then ( ^0|) implies that the decay rate is 



8 



H 2 



(l+42cg)# 



H where c. ~ — 
s 21 



(54) 



As can be seen from (^2j), the temperature behaviour is 



T 
f 



7 



H 



(55) 



After the initial stage in the nonrelativistic regime, the creation of ultrarelativistic particles slows down, while the 
nonrelativistic particles continue to decay, and we no longer expect that h is positive. Thus T < 0, although by 
(|S5|) the cooling rate is less than the non-dissipative case (q, = 0). We expect that the temperature should decrease 
monotonically after reaching its maximum T re h- By ( |5^ ) this is the case provided that 



7 -c 2 

8 s 



r 2 < i 



(56) 



We expect that ( |5q ) is easily satished in the intermediate regime, after the initial violent rate of creation of particles 
has passed, and the two-fluid mixture evolves increasingly towards 'normal' two-fluid behaviour, for which q, <C 1 
(see equation (53) of p2[ ). 



(iii) Ultra-relativistic regime 



In this last stage of reheating, the energy density becomes dominated by the radiation fluid, i.e. n 2 T n\m and 



Then 



Therefore from (Mil) we have 



implies the decay rate: 



H + 2H A 



H 



1 



(57) 



Q « 36c; 



n 2 T 



H 



The temperature change follows from Q42|) as 



T 
T 



1 - 



H where c 2 



H 



(58) 



(59) 



As argued above, we expect that Cb < c s is easily satisfied, so that the temperature continues to fall, although still 
at a reduced rate relative to the non-dissipative case. Towards the end of reheating the cosmic medium approaches a 
perfect, relativistic fluid with vanishing so that c 2 must tend to zero. For any nonzero c 2 the decay rate Q diverges 
in the limit n\m <C n 2 T. 



B. Perturbations of the expansion rate 



It is possible, given a priori the forms of the relaxation time r and the decay rate Q, to calculate the perturbations 
of the expansion rate H due to causal viscous and reaction effects, via (flO]). We illustrate this with a simple model, 
which in particular can accomodate oscillations in H during the initial stage of reheating (compare p241). 

The simple model is based on the ansatz that r _1 and Q are proportional to the expansion rate in the initial stage 
of reheating, i.e. 

t- 1 =vH , Q = (3H , (60) 

where v and (3 are positive constants (with v > 1 for a consistent hydrodynamic description). The ansatz for Q is 
consistent with (^5|) if c& is constant. Using (|6(]) in ( ff0| ) in the nonrelativistic regime, we find that the evolution of 
the expansion rate is governed by 
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H + (5 + v)HH + [3 + \v + g/3 - |cg] ff 3 « . (61) 

Now we know that 

2 2 
i? = h h where \h\ < — . (62) 

Substituting (^32|) into ( |6l| ) and linearising, we find that the perturbation h is governed by 

h + [|(5 + i/)] r 1 //. + [i (4+ 81/ + 5/3 - 36c 2 )] t~ 2 /i 

« (36c 2 - 5/3)] i~ 3 . 

This is the equation of a forced damped oscillator, as is readily seen after the change of variable to s = hit: 

h" + [§(7 + 2v)) ti + [| (4 + 81/ + 5/3 - 36cg)] /i 

« (36c 2 - 5/3)] e~ a . (63) 

The simplest case is (|45|), i.e., 5/3 = 36c 2 , for which ( |63| ) leads to the overdamped perturbation 

^(i)« £l <- 2 + £2 i-( 1+2 ^/ 3 . (64) 

Damped oscillations of h about the zero-order solution 2/(3t) occur when 

30/3 > 216c 2 + (2i/ - 5) 2 , (65) 

which shows that Cb is a damping factor, while the decay coefficient contributes to oscillation. A further constraint 
on the thermodynamic parameters arises from the requirement that the particular integral of ( |63| ) must be small 
compared to 2/(3t): 



36c 2 - 5/3 



< 1 . (66) 



5/3 - 36c 2 + 4z/ - 4 

If ( |65| ) and (|6^) are satisfied, then the damped oscillatory perturbation is given by 

h{t) pb fSt" 1 +<-( 7+2,y )/ 6 [e 1 cos(wlni) + £ 2 sin(wln<)] , (67) 

where 

= i [30/3 - 216c 2 . - (2z/ - 5) 2 ] 1/2 
is the frequency of oscillation. A simple choice that satisfies (^) and ( |66| ) is 

7/ = 5 — _2_ i 36 -2 

2 ' ^ 165 T 5 f> ' 

in which case 



rr 2 1 1 / Vll, \ 

ij rj h^Eexp ! at . 

3t 450i i 2 P ^ 33 y 

Finally, we note that the perturbations of the Hubble rate induce perturbations of the temperature via the tem- 
perature evolution equation (|42|). Using ( |60| ) and n-i 3> 7i2, Pi + P2 + tt s» niT, we find 

f + 2HT m \(3mH . 

Denoting by T the dominant zero-order part of the temperature, which is determined by H = 2/(3i), and by ST the 
perturbation induced by SH = h, this equation leads to 

(ST)' + 2HST& (±Pm-2T)h . (68) 

Equation ( |68| ) determines the temperature perturbation explicitly in terms of the Hubble rate perturbation. Clearly, 
if h has an oscillatory component, then so will ST. 
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V. CONCLUDING REMARKS 



A qualitative analysis based on a thermodynamic two-fluid model subject to causality has produced overall features 
which arc consistent with expectations. The decay rate of nonrelativistic particles implied by causality is positive 
throughout the reheating process (thus providing a consistency check on our approach), and proportional to the 
expansion. The overall behaviour of the temperature and decay rate are essentially determined by q,, the dissipative 
contribution to the sound speed. Provided that Q, exceeds the adiabatic contribution c s in the beginning of reheat- 
ing, equivalently, provided that the total number density grows by virtue of super-abundant particle creation, the 
temperature rises very rapidly at the start of reheating. It quickly reaches a maximum value (fl9|), whereafter it falls 
with expansion, at a reduced rate (determined by ct,) relative to the non-dissipative case. Large amounts of entropy 
are generated in the early stage, as shown by (Q). A simplified model ( |60| ) allows us to calculate explicitly the de- 
caying perturbations of the Hubble rate, including a case of oscillatory perturbations ( |67| ) around matter-dominated 
behavior. 

We have thus shown the capacity of a phenomenological model based on causal relativistic thermodynamics to 
predict the expected basic features of reheating with economy and simplicity, and independent of detailed knowledge 
of the interaction. This should be counted as another success of the Israel-Stewart theory (as adapted to deal with 
interacting fluids), and another argument in favour of the theory with its causality and stability properties. 
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